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ABSTRACT 


An  investigation  is  made  of  the  structure  of  the  electro- 
magnetic field  near  the  focvis  of  an  aplanatic  system  which  images  a  point 
source.  First  the  case  of  a  linearly  polarized  incident  field  is  examined 
and  expressions  are  derived  for  the  electric  and  magnetic  vectors  in  the 
image  space.  Some  general  consequences  of  the  formulas  are  then  discixssed. 
In  particular  the  synmetry  properties  of  the  field  with  respect  to  the  focal 
plane  are  noted  and  the  state  of  polarization  of  the  image  region  is  inves- 
tigated. The  distribution  of  the  time  averaged  electric  and  magnetic  energy 
densities  and  of  the  energy  flow  (Poynting  vector)  in  the  focal  plane  is 
studied  in  detail  and  the  results  are  illustrated  by  diagrams  and  in  a  tabu- 
lated form  based  on  data  obtained  by  extensive  calculations  on  an  electronic 
computer.  The  case  of  an  unpolarized  field  is  also  investigated. 

The  solution  is  not  restricted  to  systems  of  low  aperture,  and 
the  computational  results  cover,  in  fact,  selected  values  of  the  angxilar 
semi-aperture  a  on  the  image  side,  in  the  whole  range  0  <  a  <  90°.  The 
limiting  case  a  ->  0  is  examined  in  detail  and  it  is  shown  that  the  field 
is  then  completely  characterized  by  a  single,  generally  conrplex,  scalar 
function,  which  turns  out  to  be  identical  with  that  of  the  classical  scalar 
theory  of  Ali^,  Lommel  and  Struve. 

The  results  have  an  immediate  bearing  on  the  resolving  power  of 
image  forming  systems;   they  also  help  to  imderstand  the  significance  of 
the  scalar  diffraction  theory,  which  is  customarily  en^jloyed  without  proper 
Justification,  in  the  analysis  of  images  in  low  aperture,  systems. 
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1.  Introduction 

A  knowledge  of  the  structure  of  an  electrcanagnetic  field  in  the 
region  of  focus  is  of  considerable  theoretical  as  well  as  practical  interest. 
As  already  indicated  in  a  previous  report  [Wolf,1958],  information  about  the 
structxire  of  this  complex  region  is  particularly  desirable  in  connection  with 
the  design  and  the  analysis  of  performance  of  optical  systems  of  wide  angvilar 
aperture,  both  in  the  field  of  visible  and  microwave  optics.  The  construc- 
tion of  large  peo^boloids  for  use   in  reidio  astronomy  has  made  this  pzx)blem 
also  of  topical  interest. 

Apart  from  the  usual  scalar  treatments,  the  limitations  of  which 
were  mentioned  in  the  previous  report,  the  only  earlier  treatments  of  this 
problem  appear  to  be  those  of  Ignatowaky  (1919,  1920)  and  Hopkins  (19^5, 
19^5).*  The  investigations  of  Ignatowsky  have  gone  a  considerable  way  towards 
the  solution  of  the  problem.  He  obtained  formulas  for  the  electric  and 
magnetic  field  vectors  in  the  image  region  of  an  aplanatic  system**  and  also 
of  a  paraboloid  of  any  angular  aperture.  Unfortunately  his  deductions  froia 
these  formiilas  were  chiefly  confined  to  the  study  of  the  energy  flow  across 
the  central  bright  nucleais  of  the  Image;  the  electric  energy  (which  is  pre- 
simiably  what  most  detectors,  e.g.  a  photographic  plate,  record)  is  not 
discTJSsed.  The  researches  of  Hopkins  were  mainly  concerned  with  the  modifi- 
cation which  the  Airy  pattern  undergoes  aa   the  angular  apertxire  of  the  system 
is  increased.  While  his  analysis  is  not  based  on  the  full  Maxwell's  equa- 
tions, it  does  take  into  account  the  vectorieuL  nature  of  the  problem  and  his 
results  are  of  practical  interest  in  connection  with  optical  systems  the 
angular  semi -apertures  of  which  do  not  exceed  about  4o°* 

*Since  the  present  investigation  was  carried  out,  two  related  papers  appeared. 
Burtin  (1956)  determined  the  distribution  of  the  electric  energy  density  in 
the  focal  plane  of  a  system  with  angular  semi -aperture  a  =  ^5°  on  the  image 
side,  when  the  wave  entering  the  system  is  linearly  policrized.  Focke  (1957) 
considered  an  unpolarized  wave  and  studied  the  energy  density  and  the  energy 
flow  at  points  in  the  focal  plane  in  systems  of  selected  anguleur  ai)ertiires. 

**By  aplanatic  system  we  mean  one  which,  for  a  specified  axial  position  of 
the  object  point,  is  stigmatic  and  obeys  the  Abbe  sine  condition. 
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In  the  present  paper  a  thorough  investigation  is  made  of  the 
structure  of  the  electromagnetic  field  near  the  focus  of  an  aplanatic 
system  which  images  a  point  source  at  infinity.  Formulas  are  derived  for 
the  electric  and  magnetic  vectors  in  the  image  space,  and  a  number  of 
general  properties  of  the  electromagnetic  field  are   deduced,  both  for  the 
case  of  polarized,  as  well  as  impolarized  incident  waves.  The  foimulaB 
are  evaluated  for  a  large  nvmiber  of  the  basic  parameters  of  the  problem  and 
the  results  are  presented  in  the  form  of  diagrams  and  in  tabulated  form. 
Systems  of  selected  angular  semi-aperture  a  on  the  image  side  are  considered, 
up  to  the  limiting  ceise  a  ->  90  .  The  other  limiting  case,  a  ->  0,  is  exam- 
ined in  detail  and  it  is  foimd  that,  in  this  case,  the  vector  solution  is 
completely  characterized  by  one  (generally  complex)  scalar  function,  which 
is  found  to  be  identical  with  that  of  the  classical  analyses  of  Airy  (I855), 
Lommel  (1885)  and  Struve  (1886). 

The  results  have  an  immediate  bee^ring  on  the  resolving  power  of 
systems  with  high  angular  aperture;  they  also  help  our  undei^tanding  of  the 
significance  of  the  usual  scalar  methods  of  optical  diffraction  theory. 

Finally,  it  may  be  mentioned  that  our  basic  fonntilas.  Equations 
(2.30)  and  (2.51),  are  in  agreement  with  those  of  Ignatowsky  (1919);  ovir 
deductions  go  considerably  further,  however. 

Some  preliminary  results  of  this  investigation  were  reported  in 
two  previoxis  notes  JRichards  (1956a),  Richards  and  Wolf  (1956)]]. 


2.  Expressions  for  the  field  vectors  in  the  image  space  of  an  aplanatic  system. 

Consider  an  optical  system  of  revolution  which  images  a  point  sovirce. 
The  imaging  will  be  assumed  to  be  aplanatic,  i.e.  axlally  stigmatic  and  obey- 
ing the  sine  condition.  The  source  will  be  assumed  to  be  at  infinity  in  the 
direction  of  the  axis,  and   to  begin  with  it  will  be  assumed  that  the  source 
is  monochromatic,  and  that  it  gives  rise  to  a  linearly  polarized  wave  in  the 
entrance  piipll  of  the  system.  The  case  of  an  impolarized  wave  will  be  consi- 
dered later  {y    5),  by  averaging  over  an  possible  states  of  polarization. 
It  is  assumed  that  the  linear  dimensions  of  the  exit  pijpil  are  large  compared 
with  the  wavelength. 


Expressions  for  the  field  in  the  image  region  of  the  system  may- 
be derived  by  an  application  of  formulas   (2.22)   and  (2.23)    of  a  previous 
report    [Wolf,    195tiJ  •      These  formuleis  express  the  time  iudependent  parts, 
e  and  h,of  the  electric  and  magnetic  vectors 
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(2.1)  E(x,y,z,t)   =  /C{e(x,y,z,  )e'^      )   ,    H(x,y,z,t)    =  y^ -^h(x,y,z, ) 


■icotl 


at  any  point  P(x,y,z)  in  the  image  space,  which  is  not  too  close  to  the  exit 
pupil,  in  the  forn. 


(2.2)     e(x,y,z) 


-i// 


*—   exp  < ik 


^(s  ,s  )+s  x+s  v+s  z 
^^  x'  y   X   r   z 


ds  ds 
X  y 


n. 


(2.5)     h(x,y,z)  =  -  g 


ik  ((    ^^V%^     )  Vt 

^-^i-  exp  <ik  |(s^,Sy)+s^x+s^+s^z 

n 


ds  ds 
X  y 


vhere  -c  denotes  the  real  part.  Here  a  and  b  are  the  'strength  factors'  of 
the  unperturbed  electromagnetic  field  E^^^  = -^{e^^^e"""*]  ,  H^^^  =£  h^^^e"^"*^  } 
which  is  incident  on  the  exit  pupil,  i.e.  a  and  b  are  defined  by  the  relations 


(2A) 


(i) 


a  x\s.S 


yv2 


-(i)  ^  b    ^ik^ 


/R,R, 


12 


where  R     and  R     are  the  principal  radii  of  curvature  of  the  associated  wave- 
front  and  ^5'    is  the  eikonal  function.     Further 


(2.5) 


c 


2jt 
X 


is  the  vacuum  nianber,  and  X  the  vacuum  wave  length,  and  it  is  assumed  that 
the  refractive  index  of  the  image  space  is  unity.  Further  J  is  the  aberra- 
tion function  of  the  system,  s  ,  s  ,  s  are  the  components  of  the  unit 

X   y   2 


k  - 


vector  s  (vith  its  positive  direction  in  the  direction  of  propagation  of 
the  light)   along  a  typical  ray  in  the  image  space,    and  O.   is  the  solid 
angle  formed  by  all  the  geometrical  rays  which  pass  through  the  exit, 
pupil  of  the  system.     According  to  Equation  (2.2i<-)   of  the  previous  report, 
the  sti^ngth  factors  are  related  by  the   formula  (assiming  that  in  the  image 
space   t  ~  ^1  ^l) 

(2.b;  b  =   8  /s  a 

In  the  problem  under  consideration,  the  imaging  is  aplanatic  so 
that  the  wave fronts  in  the  image  space  are  spherical,  with  a  common  centre 
at  the  Gaussian  image  point.  Hence,  for  all  vectors  s  in  the  solid  angle  -H-, 


(2.7)         f   (s^,  s  )  =  0 


y 

The  strength  factors  a  and  b  may  be  determined,  with  an  accuracy 
sufficient  for  our  purposes,  by  tracing  rays  through  the  system  up  to  the 
region  of  the  exit  pupil,  and  by  making  use  of  the  laws  relating  to  the 
variation  of  the  field  vectors  along  each  ray.  Let  AQ  be  a  typical  ray 
entering  the  system  at  a  height  h  from  the  axis  and  let  0  be  the  angle 
which  the  corresponding  ray  QO  in  the  image  space  makes  with  the  axis 
(Figure  l).  Since  the  field  in  the  object  space  is  linearly  polarized, 
the  field  obtained  on  refraction  at  the  first  surface  cr  is  also  linearly 
polarized*;  and  if  the  angle  of  the  incidence  at  a:  is  small,  the  angle 
which  the  direction  of  vibrations  of  the  electric  (and  also  of  the  magnetic) 
vector  makes  with  the  meridional  plane  of  the  ray  (the  plane  containing  the 
ray  and  the  axis  of  the  system)  will  be  effectively  unchanged  by  refraction.** 


*     The  direction  of  polarization  will, of  course, be,  in  general,  different 
for  different  rays. 

**  The  truth  of  these  statements  follows  from  Fresnel  formulas  for  refraction 
of  a  plane  wave  on  a  plane  interface  between  two  homogeneous  medis  of  dif- 
ferent refractive  indices.  Within  the  accuracy  of  the  present  approxima- 
tion (geometric  optics  limit,  i.e.  k  ->  oo),  these  formulas  also  apply  to 
the  present  case  \rtien  neither  the  wavefronts  nore  the  refracting  bound- 
earies  are  necessarily  plane. 
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FIGURE    I 


For  caption  to  Figure  see  page  1^0, 
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Moreover  [see  Bom  and  Wolf  (1959) >  S3«l«5]>  in  an  homogeneous  medium  the 
direction  of  vibration  remains  unchanged  along  each  ray.  By  successive 
repetition  of  these  arguments  from  surface  to  surface  it  follows  that  the 
field  in  the  region  of  the  exit  pupil  is  also  linearly  polarized^,  and  that, 
provided  the  angle  of  incidence  at  each  surface  of  the  system  is  small,  the 
angles  which  the  e  and  h  vectors  in  the  exit  pupil  make  with  the  meridional 
plane  of  the  ray  are  equal  to  the  corresponding  angle  in  the  object  space. 
Further,  the  magnitude  of  the  field  vectors  in  the  region  of  the  exit  pupil 
may  be  determined  frcm  the  gecanetrical  law  of  conservation  of  energy,  taking 
into  account  the  fact  that  the  system  obeys  the  sine  condition. 

Let  iT  be  the  'focal  sphere',  i.e.  the  sphere  with  centre  at  0 
and  with  radius  f  equal  to  the  focal  length  of  the  system.  Then,  accord- 
ing to  the  sine  condition 

(2.8)  h  =  f  sin  0  . 


This  relation  implies  that  the  emergent  ray  meets  the  focal  sphere  at  the 
same  height  at  which  the  corresponding  ray  in  the  object  space  entered  the 
system  (see  Figure  l). 

Consider  now  all  the  rays  which  enter  the  annulus  bounded  by 
circles  of  radius  h  and  h  +  Sh.  Let  6S  be  the  area  of  the  annulus  and 
6S,  be  the  corresponding  area  of  the  focal  sphere,  and  let  e  and  e  be  the 
electirLc  amplitude  vectors  on  SS  and  6S^  respectively.  Then  we  may  write 


o    -    ^  -  g    1 


(2.9)         %-KK^         °  '   ^1  =  ^o  ^1 


> 


where  e  and  e.  are  unit  vectors,  in  the  direction  of  e  and  e.,  respectively, 
o      1  o      1 

/  and  /  are  amplitude  factors  and  <^  and  S\    are  the  corresponding  values 
of  the  eikonal.  Since  the  field  is  linearly  polarized,  /  ,/. ,  e  and  e  may 
be  taken  to  be  real,  provided  that  the  origin  of  the  phase  is  suitably  chosen. 
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According  to  the  geometrical  optics  intensity  law, 

(2.10)  ^^8S   =-^?6S,  , 
^  0   0      11' 

where  we  assumed  that  the  refractive  index  of  the  object  space,  like  that 
of  the  image  space  is  unity,  and  the  losses  of  energy  due  to  reflection 
and  absorption  within  the  system  are  negligible.  Now  from  the  figure, 

(2.11)  5S   =  5S,cosO  , 
^  o      1     ' 

so  that  according  to  (2.10) 

(2.12)  ^^  =  -/  /cosG  . 

Now  e    may  be  identified  with  the  vector  e  ^-^^   of  (2.4),  with  R-    =  R     =  f 
and  -^    =  ■^^    ,  Hence 


(2.15)  a     =  t£^JcosQ  e 


To  determine  e,  it  is  convenient  to  introduce  two  unit  vectors 

g  and  g^  in  the  meridional  plane  of  the  iray,  such  that  g  is  perpendicu- 

lar  to  the  ray  in  the  object  space  and  g^  is  perpendicvilar  to  the  ray  in 

the  image  space,  and  both  are  directed  away  from  the  axis  (see  Figure  1). 

Let  Ox,  Oy,  Oz  be  Cartesian  rectangtilar  coordinate  axes,  with  origin  at 

the  Gaussian  focios,  with  Ox  in  the  direction  of  e  and  with  Oz  along  the 

o 

axis  of  the  system,  pointing  away  from  the  plane  of  the  exit  pupil  into 
the  image  space.  Finally,  let  i,  j  and  k  be  unit  vectors  in  the  direction 
of  the  coordinate  axes. 

The  electric  and  the  magnetic  vectors  are  orthogonal  to  the  ray 
[see  Bom  and  Wolf  (l959),§  5«l]«  Hence  e  lies  in  the  plane  of  g,  and 
g^^  s  ,  i.e. 


A 


(2.11+)        e^  =  og^  +  p(  g^^  s  ) 


where  a  and  p  are  some  constants.  To  determine  a  and  p  ve  make  use  of 
a  resxilt  mentioned  eeirlier,  namely  that  as  the  light  traverses  the  system, 
the  angle  between  the  electric  (and  also  the  magnetic)  vector  eind  the  meri- 
dional plane  of  the  ray  remains  constant.  Moreover,  each  of  these  vectors 
remains  on  the  same  side  of  the  meridional  plane.  These  results  imply  that 


Si  •  ®1  =  s^  •  ^ 


(2.15)   < 


/^ 


(g^A  s^)  .  e^  =  (g^A  k)  .  i  . 
On  substituting  from  (2.1^)  into  (2.15)  we  find  that 
(2.16)    a  =  g^  .  i  ,  P  =  (g„Ak)  .  i  =  g^  .  (kAi)  =  L  .  j  , 


5^.1,3  =  (g^Ak)  .  i  =g^ 
and  from  (2.12),  (2.15)  and  (2.15)  it  follows  that 


(2.17)    a  =  f,i^^yHHio[(g^.  i)|^  +  (g^.  j)(g^At)]  . 

It  will  now  be  convenient  to  introduce  spherical  polar  coordinates 
r,  ©,  0  (r  >  0,  0  <  e  <  jt,  0  <  0  <  2k),   with  the  polar  axis  0  =  0  in  the 
z-direction,  and  with  the  azimuth  0=0  containing  the  electric  vector  in  the 
object  space.  The  ccanponents  s  ,  s  ,  s  of  the  imit  vector  s  along  a  ray  in 
the  image  space  and  the  coordinates  (x,  y,  z)  of  a  point  P  in  the  image 
region  may  be  expressed  in  the  form 

(2.lBa)   a  =  sin  Q  cos  0  ,   s  =  sin  0  sin  0  ,   s  =  cos  9  , 

(2.18b)   ^  =  Tp  sin  Op  cos  0p  ,  y  =  rp  sin  9p  sin  0p  ,  z  =  Tp  cos  ©p  > 

so  that  the  term  in  the  exponent  of  the  Integrals  (2.1)  and  (2.5)  becomes 

(2.19)  s  X  +  s  y  +  8  z  =  Tp  cos  €  , 
where 

(2.20)  cos  e  =  cos  0  cos  Op  +  sin  0  sin  0  cos  (0  -  0  )  . 
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Let  (0,0)  and  ( 0. ,  0^ )  be  the  polar  angles  of  g  and  g^ 

respectively.  Then  evidently  (see  Figure  l) 

(2.21)    0^  =  0^  =  0  -  rt  ,  0^  =  I  ,  0^  =  I  .  ©  , 


so  that 


(2.22) 


5^  =  -  COS0  i  -  sin0  j  , 

5^  =  -  COS0COS0  i  -  co60sin0  j  +  sin©  k  • 


It  follovs  on  substitution  from  (2.22)  and  (2.l8a)  into  (2.1?)  and  (2.6) 
that  the  Cartesian  components  of  the  strength  vectors  a  and  b  are 

a  =  tX    ycosOfcosO  +  sin  0(1  -  cos0)l  , 

jC  O  — ' 

(2.25)        ^     %  =  f^/cos©[(cos©  -  I)co608in0]    , 
a     =~t^  ycos0sin©cos0  , 

b     =  f  j2  /cos©  [( cos©  -  I)cos08in0]   , 
(2.2if)        \   \  =  ^J^  /co8©[l  -  8in^0(l  -  cos©)^)    , 


b     =-tJ^  7cos0sin0sin0  . 
V^    s  o 


Finally,  we  also  need  an  expression,   in  terms  of  0  and  0,   for 

the  quantity  ds     ds  /s   .  which  enters  our  basic  diffraction  integrals 
X       y     z 

(2.2)  and  (2.5).     According  to  equation  (2.6)  of  the  previous  report,   this 
quantity  represents  the  element  diT.  of  the  solid  angle  and  is  given  by 


ds  ds 
(2.25)  ^     y     =     d-O-     =  sin0d0d0  . 
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On  substituting  from  (2.25),  (2.2^),  (2.25),  (2.11)  and  (2.7) 
into  (2.2)  and  (2.5)  we  obtain  the  following  expressions  for  the  Cartesian 
components  e  and  h  : 


e  = 


a  2n 
—  J    I  VcosOsin©  <cos»+fl-co80)B-"-  " 
o  o 


y cos 08 in©  \  cos0+(l-co80)Bin  0/  exp[ikr  cose  Jd0d0  , 


a  2jt 


(2.26)  ■/     e  =       —       I  ycos©sine(l-co80)cos0sin0exp[ikr  cose  Jd9d0  , 


o  o 


a  2jt 


e  =       —   I    /  y  cos  ©sin  ©co80exp  fikTpCose  Jd©d0  ; 
o  o 


a  2jt 


h  = 

X 


lA 


(2.27) < 


—  ycos©sin©(l-cos©)cos0sin0exp[ikr  cose  J  d©d0  , 

o  o 

a  2jt 
iA    r  r  r—r:.^   J,    /,         _x    .    2, 


h  = 7)1  ycosGsin©  <l-(l-cos©)8in  0V  exp[ikr  cose  ]d©d0  , 


h  = 

z 


o  o 

a  23t 
±A    f  r  /■rrr^_._2 
rt 
o  o 


—       I  /coiesin  ©sin0exp(ikrpcose  2d©d0  . 


Here  cos  e   is  given  by  (2.19),  a,  is  the  euagiilar  semi-aperture  on  the  image 
side,    i.e.   ax  is  the  angle  which  the  diameter  of  the  exit  pupil  subtends  at 
the  gecanetrical  focus  «ind  A  is  the  constant 


(2.28)         A     =  kf>^^  /  2     =     %tj^^  I  X  . 
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The  integration  with  respect  to  0  can  irmnediately  be  carried 
out  vith  the  help  of  the  follovfing  formulas*  vhich  are  valid  for  any 
integral  value  of  n  : 


2rt 

'  cos  n  0  exp[i  p  cos   (0  -   7)]'i0  =  2«  i^  J^^  (p)    cos  n  7  , 


J 

o 


(2.29)    < 


2rt 


/  X\ 

sin  n  0  exp[i   p  cos   (0  -  7)]d0  =  2rt  i     J^^  (p)   sin  n  7 


o 


Here  J  (p)  is  the  Bessel  function  of  the  first  kind  and  order  n.   If  in 
(2.26)  and  (2.27)  we  use  the  identities  cos  0  sin  0  =  -g  sin  2  J?, 
8in^0  =  i  (1  -  cos  20)  and  apply  (2.29)  we  finally  obtain  the  following 
expressions  for  the  components  of  the  field  vectors  at  a  point  P  in  the 
image  jregion: 


These  formulas  may  be  derived  as  follows:  We  start  from  the  integral 


representation  of  J 


2n 


I  ei(°i^  + 


pcos0) 


d0  =  ati^J  (p) 


(of.  Watson  (1952),  p.  20  (5)).  We  change  0  into  0-7  and  multiply 
both  sides  of  (2.5O)  by  e  ''  and  express  the  resulting  formula  as  follows,* 


2s.  2rt 

pcos(0-7)^^  ^  ,    (  «._«uipco3(0-7) 


I  cosn0e 
o 


'd0  +  i  I  8inn0e^ 
o 


d0»  2ni  J  ( p) [cosn7  +  isinnTJ  . 


Each  side  consists  of  two  tenns,  one  of  which  is  an  even  function  of  7  and 
the  other  an  odd  function  of  7.  This  is  only  possible  if  the  even  terns 
are  equal  to  each  other  and  so  are  the  odd  terms,  and  this  in5>lies  (2.29). 
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e^(P)  =  -  iA  (I^  +  Ig  cos  20p)  , 


(2. JO)  <   e^CP)  =  -  iA  I2  sin  20^  , 


e  (?)  =  -  2A  I  cos  0  , 

Z  X  r 


h^(P)  =  -  iA  I2  sin  20p  , 


(2.51)  <  hy(P)  =  -  iA  (I^  -  I2  cos  20p)  , 


h^(P)  =  -  2A  I^  sin  0p  , 


where 


1=1  (kr  ,©  ,a)  =  I  \fcoB~9   sin  6  (1  +  cos  0)  J  (kr  sin  0  sin  0  ) 

o 

.    exp    (jLkTp  cos  0  cos  0^   d0  , 


a 


(2.52)     < 


I^  =  I^(krp,Op,a)   =  I  v/cos  0  sin^0  J-^C^rp  sin  0  sin  0p) 

o 

•  exp    [ifcTp  cos  0  cos  oJ]    d0  , 


a 


I^  =  l2(krp,0p,a)   =  j  v/cos  0  sin  0  (l  -  cos  ©)  J2(^rp  sin  0  sin  0p) 

o 

♦  exp  [lltTp  cos  0  cos  0p]  d0  . 


The  formulas  (2.3O)  and  (2. 51)  represent  the  analytic  solution 
of  our  problem.  They  express  the  field  at  any  point  P  (spherical  polar 
coordinates  r  ,  0  ,  0  )  in  terms  of  the  three  integrals  I  ,  I^  and  I  . 
We  shall  now  discuss  some  consequences  of  these  formulas. 
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3-   The  Image  field 

It  is  convenient  at  this  stage  to  introduce  certain  'optical 
coordinates',  which  are  a  natural  generalization  of  the  coordinates  (defined 
by  (3 -la)  below)  employed  frequently  in  connection  with  diffraction  in  systems 
with  low  angular  aperture.  We  define  these  optical  coordinates  by  the  formulas* 


(5.1)  < 


2         2 
u  =  kr_  cos  0  sin  a  =  kz  sin  a  , 


/~2    2~ 
V  =  krp  sin  0  sin  a  =  kv/x  +  y   sin  a 


From  now  on  we  shall  omit  the  subscript  P  in  the  symbol  0  for  the  azimuthal 
angle,  and  specify  the  point  P  of  observation  by  the  three  parameters,  u,  v 
and  0(  u—  0,  v>0,  O<0<2rt).  The  geometrical  focal,  plane  is  given  by 
u  =  0,  the  ajcis  by  v  =  0  and  the  edge  of  the  geometrical  shadow 

L/x  +  y  =  +  z  tan  a  j  by  v  =  +  u  sec  a  . 

The  integrals   (2.52)   are  now  regarded  as  functions  of  u  and  v. 


I„(u,v) 


/ 


v^cos  0  sin  0  (1  +  cos  0)  J  (^„f^°^  |exp 


o\   sin  a 


a. 
(3.2)  J        I^(u,v)  =|y^5r0  6in20  J^t^fiS-^jexp 


a 


l2(u,v) 


V 


iu  cos  0 
sin  a   -^ 


V  sin  0 


ycos  0  sin  0  (1  -  cos  ©)  J  ( — jexp 


iu  cos  0 

^  2    , 
sin  a   -■ 


d0  , 


iu  cos  0 

4  2    _, 

1—  sin  a  -* 


d0  , 


d0 


^When  the  angular  aperture  is  small  (a«  l),  the  formulas  (5'1)  reduce  to 


(5.1a) 


u  .^  k  (  a/R  )  z  , 


k  (  a/R  )  s/x  +  y^  , 


where  sin  a  /^  tan  a  =  a/R,  a  being  the  radius  of  the  exit  pupil  and  R  the  dis- 
tsuice  between  the  exit  pupil  and  the  focal  plane. 
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We  note  that 


(5.3)     Iq  (-u,v)  =  I^  (u,v)   ,   (n  =  0,1,2)  , 


where  the  asterisk  denotes  the  complex  conjugate.  From  {2.^0),    (2.3I) 
and  (3.5)  we  note  the  folloving  relations  which  exist  between  the  com- 
ponents of  the  field  vectors  at  any  two  points  P^(u,v,0)  and  P^(-u,v,0), 
vhich  are  symmetrically  placed  with  respect  to  the  focal  plane: 


e^  (-u,v,0>  =  -  e^  (u,v,0)  ,   h^  (-u,v,0) 


h^  (u,v,0)  , 


(5.^)  <   e  (-u,v,0)  =  -  e*  (u,v,0)  ,   h  (-u,v,0)  =  -  h  (u,v,0)  , 


e^  (-u,v,0)  =  e^  (vi,v,0)  ,   h^  (-u,v,0)  =   h^  (u,v,0)  , 


If  |e  I  denotes  the  amplitude  and  $   the  phase  of  e  ,  the  first  rela- 
X  ^  X  X 

tion  in  (5.I+)  implies  that 

(5.5a)    I  e^  (-u,v,0)  1  =  I  e^  (u,v,0)  |   , 
and 

(5.5b)   $  ^  (-u,v,0)=-^   (u,v,0)  +  jt    (mod  2n)   . 


The  relation  (5.5a)  shows  that  for  any  tvo  points  which  are  symmetrically 
situated  vith  respect  to  the  focal  plane,  the  aniplitude  of  the  e  's  are 


*The  quantity  laod  2n  on  the  right  of  an  equation  denotes  that  the  two 
sides  of  the  eqxiation  are  indeterminate  to  the  extent  of  an  additive 


constant  2mjr  vhere  m  is  any  integer. 
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the  same,  and  (5.5b)  shovs  that  there  exists  eJ.8o  a  slisple  relation 

between  the  phases.  The  appearance  of  the  additive  factor  it  on  the 

right  of  (5.5b)  is  connected  with  the  veil-known  phase  ancHnaly  at 

focus  [cf  Bom  and  Wolf,  (1959)  ,  S   8-8>  e^l-  (^5)]  •  From  {^.h)   we 

see  that  relations  of  the  form  (5.5a)  holds  also  for  all  the  other 

components.  However,  relation  of  the  form  (5.5b)  holds  for  the  phases 

of  e  ,  e  ,  h  and  h  only.  For  the  phase  of  the  remaining  two  com- 
X    y    X     y 

ponents  we  have,  in  place  of  (3.5b),  relations  of  the  following  form: 
(5.5c)    $^  (-u,v,0)  =-$'^  (u,v,0)    (mod  at). 

Thus  the  components  of  e  and  h  in  the  direction  of  the  axis  of  revo- 
Itition  of  the  system  have  no  phase  anomaly. 

From  (2.50)  and  (2.51)  it  is  also  seen  that 


\   (u,v,0) 


(5.6)   <    hy  (u,v,0) 


\   (u,v,0) 


ey  (u,v,0  -  |)  , 


e^  (u,v,0  -  |)  , 


e^  (u,v,0  -  |)  . 


Hence  in  any  fixed  plane  of  observation  (  u  =  constant  ),  the  e  and  h 
fields  are  the  same  but  are  rotated  with  respect  to  each  other  by  90 
around  the  z-axis;  this,  of  coxirse,  might  have  been  expected,  since  the 
e  and  h  fields  in  the  object  space  have  this  relationship,  and  the  laws 
relating  to  the  transmission  of  these  fields  through  the  system  are  the 
same. 


16  - 


5.1  Polarization  of  the  image  field 

To  examine  the  state  of  polarization  of  the  image  field, 
ve  separate  the  real  and   imaginary  parts  of  the  integrals  (5' 2)  emd 
write 

(5.7)     I  (u,v)  =  I  ^""^  (u,v)  +  il  ^^^  (u,v)  ,   (n  =  0,1,2) 


where  I  ^   and  I     are  real.  We  also  write 
n        n 


(5.8)     e  (u,v,0)  =  p  (u,v,0)  +  iq  (u,v,0)  , 

where  p  and  q  ajre  real  vectors;  they  are  a  pair  of  conjugate  semi- 
diameters  of  the  polarization  ellipse  of  the  electric  vector.  According 
to  (2.50),  (5.7)  and  (5.8) 


p^(u,v,0)  =  A(I^^^^+  I2^^^CO820)  ,   q^(u,v,0)  =-A(l^^^^l2^''^os20), 


(5.9)  { 


Py(u,v,0)  =  A  Ig^sin  20  , 
P  (u,v,0)  =  -  2A  1.^^'    cos  0 


qy(u,v,0)  =  -A  Ig^""^  sin  20  , 
q^(u,v,0)  =  -2A  I^^^^cos  0  . 


Since,  according  to  (5.^)  the  integrals  I^,  I^  and  I ^  are  all  real 
when  u  =  0  (focal  plane),  it  follows  that 


Pj^  (o,v,0)  =  0  , 


q^  (o,v,0)  =  -A[l^(o,v)  +  l2(o,v)  cos  20]  , 
q^  (o,v,0)  =  -A  1^(0, v)  sin  20  , 
p  (o,v,0)  =  -  2A  Ij^  (o,v)  cos  0  ,   q^  (o,v,0)  =  0  . 


(5.10)  <;  Py  (o,v,0)  =  0  , 
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From  (5.10)  we  see  that  in  the  focal  plane  p.q  =  0,  i.e.  the  conjiigate  semi- 
diameters  are  at  right  angles  to  each  other;  hence  in  the  focal  plane,  p  and 
q  are  the  semi-axes  of  the  polarization  ellipse  of  the  electric  vector.  More- 
over,  the  p-axLs  is  perpendicular  to  the  focal  plane  and  the  q-axis  lies  in 
the  focal  plane.  Thus,  the  polarization  ellipse  of  the  electric  vector  at  any 
point  in  the  focal  plane  is  at  right  angles  to  the  focal  plane.  The  angle 
'X(v,0)(-rt/2  <  "X  <  rt/2)  between  the  plane  of  the  polarization  ellipse  and  the 
plane  0=0  (the  xz-plane)  is  given  by 


q^(o,v,0)         Ip(o,v)  sin  20 
(3.11)    tanX(v,o)  =  q^(o,v,0)  =  I^(o,v)  H-  I^(o,v)  cos  20  ' 


and  the  two  axes  of  the  polarization  ellipse  are  in  the  ratio 


(5.12)    p  (v,0)  = 


Ip^  (o7V,0)| 


J^o,v,0)   +  q^(o,v,0) 


2|l^(o,v)  cos  01 


y  I?(o>v)  +  1^(0, v)  +  21  (o,v)l  (o,v)cos20 


Along  the  y-axis  (0  =  ^  or  ^)  only  one  of  the  components  in  (3. 10), 
namely  a,   is  different  from  zero.  Hence  at  each  point  of  the  y-axis  and,  in 
particular,  at  the  focvis  itself,  the  electric  field  is  linearly  polarized  in 
the  x-direction,  i.e.,  in  the  direction  of  vibrations  of  the  electric  field 
in  the  object  space.  Since  p  also  vanishes  ^en  I  (o,v)  =  0,  the  electric 
field  is  also  linearly  polarized  along  circles  centered  on  the  focus,  whose 
radii  (measured  in  v-units)  are  given  by  the  roots  of  this  equation;  however, 
along  these  circles,  the  direction  of  vibrations  is  not  in  general  in  the 
x-direction. 
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The  behavior  of  X  and  p  in  the  focal  plane  of  an  aplanatic  system 
of  angular  semi-aperture  a  =  6o  is  shown  in  Figiire  2.  The  figure  vas 
computed  from  the  formiilas  (j.ll),  (5»12)  and  (5«2)  by  means  of  the  Manchester 
University  Electronic  Coiiiputor  Mark  I.  The  details  of  these  and  other  com- 
putations relating  to  this  paper  ere   given  in  a  thesis  by  one  of  the  authors 
(Richards^  1956b). 

Figure  2(a)  shows  that  in  the  central  nucleus  of  the  image  %  -^0  , 
i.e.  in  this  region  the  electric  vector  vibrates  in  planes  nearly  parallel 
to  the  plane  which  contains  the  direction  of  the  electric  vibrations  in  the 
object  space  and  the  axis  of  the  system.  According  to  Figure  2(b)  the 
axial  component  (e  )  is  not  negligible  in  all  perts  of  the  central  nucleus. 

The  state  of  polarization  of  the  field  at  points  on  the  axis  of  the 
system  (v  =  O)  is  also  of  interest.  Since  I  (u, O)  =  I  (u, O)  =  0  for  all 
values  of  u,  it  follows  from  (2.50)  that  on  the  axis  e  =  e  =0.  Hence 
the  electric  vector  at  each  point  on  the  axle   of  revolution  in  the  imflge  space 
is  linearly  polarized,  and  its  direction  is  the  same  as  the  direction  of  the 
electric  vector  in  the  object  space. 


5.2  The  energy  density 


Let 


(5.13)  < 


<vju,v,0)>  =^<t2>=  i^Ce.t*), 


<v^(u,v,0)> 

<V  (U,V,0)> 


<;Wg(u,v,0)>  +  /wju,v,0)> 
j^  (  e  .  e  +  h  .  h  )  , 
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FIGURE  2(a) 


For  captions  to  Figxires  see  page  U0« 
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FIGURE  2(b) 
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denote  the  electric  energy  density,  the  magnetic  energy  density  and  the  total 
energy  density,  each  averaged  over  a  time  interval  which  is  large  compared 
to  the  basic  period  T  =  2rt/m.  From  (2. 50)  and  (5. 15),  ve  have 


<v^(u,v,0)>  =  3^<^llj2+  k\l^\^cosH  +  \l^f^  2co820/^(l^l/)y, 


2_._2, 


^5-1^)    \     <yj^^'^'<l>)>  =  ^\K\^  MlirBin^0  +   llsi    -  2cos20/^(l^l2   ) 


<v(u,v,0)>=^^|lj2^2|lj2^   |I. 


where  ^   denotes  the  real  part. 

We  note  that  because  of  (5.5)  jas,  of  course,  is  also  evident  fron 
(5.^)11 

(5.15)  <w^  (-u,v,0)>  =  <w^  (u,v,0)>   , 

with  similar  relations  involving  <(  v  ^  and  <  w  >  .  Hence  the  distributions 
of  the  time  averaged  electric  energy  density,  magnetic  energy  density  and 
total  energy  density  are  symmetrical  about  the  focal  plane  u  =  0» 
Further  we  see  that 

(5.16)  <w^  (u,v,0)>  .  <w^  (u,v,0  -  |)>  . 


Thus  the  distribution  of  the  time  averaged  magnetic  energy  density  is 
identical  with  the  distribution  of  the  time  averaged  electric  energy  density, 
but  the  distributions  are  rotated  with  respect  to  each  other  by  90°  about  the 
axis  of  the  system.  We  also  note  that  <^w  ^  is  independent  of  0,  so  that 
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the  loci  of  constant  time  avei-aged  total  energy  density  are  surfaces  of 
revolution  about  the  axis  of  revolution  of  the  system. 

Since  for  all  u  l,(u,0)  =  I  (u,0)  =  0,  ve  see  frtan  (3.1^)  that 
along  the  axis  of  revolution  In  the  image  space* 

(5.17)    <v^>  =  <V  =  I  <">  =  I^  Ilo(^>0)l'- 
When  u  =  V  =  0,  we  have  from  (3.2) 


(5.18)    I  (0,0)  =  I  >/cos  e  sine  (1  +  cos  O)  dO 


cos  ^5/^^a  (1+1  cos  a)  >, 


so  that  at  the  focus  itself 

(5.19)  <v^>    -    <v^>    =|<v>    =33^<^    1-  i  cos^5/2)^(i  ^1  ^^3^) 

In  Figvire  3  contovirs,  computed  from  (3'!^),   of  the  time  averaged 
electric  energy  density  in  the  focal  planes  of  aplanatic  systems  of  selected 
angular  semi-ai>erture8  a  are  shown.  The  contours  of  the  time  averaged 
magnetic  energy  density  are,  according  to  (3.l6),  Identical  with  those  for 
the  time  averaged  electric  energy,  but  are  rotated  by  90  about  the  normal 
to  the  plane  of  the  figure. 


Curves  which  illustrate  the  behavior  of  these  quantities  along  the  suds  of 
revolution  are  given  in  Richards  (1956a),  358.  The  paper  aLlso  contains 
curves  which  give  the  ratio  (as  function  of  a)  between  the  values  at  focvis 
computed  from  (3*19)  and  those  computed  by  the  application  of  the  Huygens- 
Kirchhoff  scalar  diffraction  theory. 
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FIGURE    3(a) 


For  captions  to  PI  glares  see  page  IjO. 
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FIGURE   3(b) 
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FIGURE    3(c) 
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FIGURE    3(e) 
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The  first  diagram  in  Figure  5  represents  the  limiting  case  a  -^  0  . 
We  shall  see  later  (§  k)   that  in  this  case  our  solution  reduces  to  that 
obtained  on  the  basis  of  the  usual  scalar  diffraction  theory.  Thus  ^(a)  is 
the  ordinary  Airy  diffraction  pattern;  the  contours  are  circles,  the  first 
zero  contour  being  given  by  v  =  1.22rt  ^   ^.3^.     As  a   increases,  the  pattern  is 
seen  to  lose  its  rotationsil  symmetry;  the  contours  in  the  neighbourhood  of 
the  focus  axe  then  approximately  elliptical,  with  their  major  axes  in  the 
direction  of  the  electric  vector  in  the  object  space  (0=0).  Further 
away  from  the  center  the  contours  are  of  a  more  complicated  form  and  closed 
loops  appear  around  certain  points  in  the  azimuths  0=0,  n/2,   n,   5rt/2. 
For  large  values  of  a,  the  ellipticity  of  the  contours  beccmes  quite  pro- 
nounced and  the  first  minimum  in  the  azimuth  0  =  0  is  well  outside  the 
first  zero  of  the  Airy  pattern,  while  in  the  azimuth  0  =  n/2  it  is  well 
inside.  Hence,  if  the  wave  in  the  object  space  is  linearly  polarized  and 
detectors  of  electric  energy  (e.g.  a  photographic  plate)  are  used,  our  solu- 
tion predicts  an  increase  in  the  resolving  power  in  wide  apertiire  systems 
for  measurements  in  the  azimuth  at  right  angles  to  that  of  the  electric 
vector  of  the  incident  wave.  This  conclusion  is  in  agreement  with  a 
prediction  of  Hopkins  (19^5>  19^5)  and  appeeirs  to  be  supported  by  early 
experiments  recorded  by  Carpenter  (1901). 

The  behavior  of  the  time  averaged  electric  energy  density  along 
the  azimuths  0  =  0  and  0  =  jt/2  in  the  focal  plane  is  shown  in  detail  in 
Figure  k.     It  is  seen  that  along  the  azimuth  0=0,  the  minima  are  not 
true  zeros  J  the  minima  are,  however,  zeros  along  the  azimuth  0  =  jt/2. 
The  behavior  of  the  time  averaged  total  energy  density,  which,  as  already 
noted,  is  independent  of  the  azimuth,  is  also  shown  in  Figures  k. 


5.5  The  energy  flow  (Poynting  vector) 

THiS   time  averaged  Poynting  vector  is  given  by 


(5.20)        <S>  =  E7<EaH>  =  -S-^it/^h*) 
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On  substituting  from  (2.50)  and  (2.51)   into  (3.20),  we  obtain  the  follow- 
ing expressions  for  the  ccanponents  of   <  S  )    : 


<\>  = 


^cos 


0^ 


*    *. 


^1  f^2  -  ^c,  ' 


I 


cA 


(J.21)  y  <s  >  =  ^ 


Sin  i,  4!|  I^  (l/.  l/)  i  , 


<^>-  ^{  lU'  -  Uc'^ 


^ 


where  'a'  denotes  the  imaginary  part. 

We  see  that  the  magnitude  of  the  time  averaged  Poynting  vector 
is  given  by 


(3.22)     |<l>|   =  g!|(|lj2-  |l2l^)%M^[l,(l/-l/)]' 


(1/2) 


Since  this  expression  is  independent  of  the  azimuth  0,  the  loci  of  constant 
I  <C  S  ) I  are  surfaces  of  revolution  about  the  axis  of  the  system. 

Along  the  axis  of  revolTition  of  the  system  in  the  image  space 
(v  =  0)  (5.22)  reduces  to 


cA 


(3.25)    |<S  >|  =  g-  I  I^  (u,o)  I  ; 
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On  comparison  with  (5«17)  ve  see  that  the  relation  |<  S  >|  =  c<w>  holds 
at  all  points  on  the  axis. 

From  (5.21)  we  readily  deduce  that  at  each  point  in  the  image 
space  the  Poynting  vector  lies  in  the  meridional  plane  which  contains  that 
point .  and  that  it  makes  an  angle  7  with  the  positive  z-axls,  given  by 


Mh   <^2  -  ^1  ' 


Sin  7  = 


( 


11^-  \i/)^  ^k-ij^ 


h  ^h'-   ^o*> 


t 


(1/2) 


(3.21j)  S 


III"  -  Hoi" 


cos  7  = 


V 


(1^0 1 


'  - 1^/;  ^ 


m4i<^2-v] 


2l(l/2) 


Using  the  relations  (5.5)  we  see  from  (5.22)  that 


(3.25)      KS  (-u,v,0)>|  =  |<S  (u,v,0)>  I  , 


i.e.  the  magnitude  of  the  time  averaj^ed  Poynting  vector  is  symmetric  with 
respect  to  the  focal  plane.  Further,  from  (3.2^),  we  have  sin  /7(-u,v,0)j  = 
-  sin  /7(u,v,0)}  ,  cos  |7(-u,v,0)|  =  cos  |7(u,v,0)|  ,  so  that,  for  n  ^  0, 


(3.26)      7  (-u,v,0)   =  2jt  -  7  (u,v,0)  . 
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When  u  =  0,  ve  have  from  (5.21),  since  the  integrals  I^(o,v) 


(n  =  0,  1,  2)  are  all  real. 


(5.27)    <:s^>   =  <Sy>  =   0 


i.e.  at  any  point  in  the  focal  plane,  the  time  averaged  energy  flov,  &8 
represented  by  the  Poynting  vector  is  perpend! cvilar  to  the  focal  plane. 
Moreover, 


(5.28)   < 


<(S^>  >  0     if   |l^  (o,v,)|  >\I^   (o,v)|  , 

and 

<;s^>  <  0    if   |l^  (o,v  )|  <|i2  (o,v)l  . 


The  second  line  implies  that  at  points  in  the  focal  plane  where 

1 1  (0,v)|  <  1 1  (0.v)|,  the  energy  flov  is  directed  back  towards  the 

object  space.* 

The  computed  ciirves,  which  show  the  behavior  of  <(S  >  in  the 
focal  plane  of  aplanatic  systems  of  selected  euigiilar  semi-aperture  a 
are  shown  in  Figure  5.  It  is  of  interest  to  note  that  these  (normalized) 
curves  do  net  substantially  depend  on  a.  In  particular,  in  the  region  of 
the  Rayleigh  limit,  where  the  durves  do  not  drop  below  0.8  of  their  respec- 
tive maxima,  the  differences  (for  any  particular  v- value  in  this  region) 
are  less  than  1  "^/j ,  in  the  full  range  (O  <  a  <  90°)  of  the  angtiar 
semi-aperture. 


*This  interesting  fact  was  already  deduced  by  Ignatowsky  (1919). 
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1*^.  Approximate  form  of  the  solution  for  syatema  vlth  a  small  angular 
aperture. 

The  structvire  of  images  in  systems  vlth  small  ajigular  aperture 
is  tusually  investigated  on  the  basis  of  a  scalar  diffraction  theory.  The 
Intensity  distribution  in  the  focal  plane  of  an  aberration-free  system  of 
revolution  which  images  a  point  soiorce  was  first  determined  in  this  way  by 
Airy  (1855).  His  analysis  was  extended  by  Lommel  (I885)  and  Struve  (I886), 
who  determined  the  distribution  in  the  whole  three-dimensional  neighbor- 
hood of  the  focus.   It  is  of  interest  to  examine  the  form  which  our  solution 
taJces  when  the  angular  semi -aperture  a  is  small,  and  to  compare  it  with 
these  early  solutions. 

When  a  is  smeLLl  enotigh  and  u  and  v  are  not  large  con^jared  with 
vinity,  we  may  replace  the  trigonometric  factors  in  the  amplitudes  of  the 
integrands  in  I  ,  I  and  I  in  (5.2)  by  the  leading  term  in  their  series 
expansions.  In  the  exponents  of  the  integrands  we  retain  both  the  first  and 
the  second  term  in  the  expansion  of  cos  d  since,  when  u  is  not  small  com- 
pared to  unity,  the  contribution  of  the  second  term  is  evidently  also 
significant.  The  integrals  then  become 


C^-l)  < 


Iju,v) 


=  2  exp 


r  T  a 

iu 

2 
a 

u  J  o 


0  J  (— )  exp 

o  ^  a    " 


-iuO 


r  1  ct 


I-j^(u,v)  =  exp 


iu 


a 
u  -J  o 


g"  Jl  (^)  exp 


ax 


■iuO 


dO  , 


ax 


de 


r  1  a, 


\{^,^) 


2  exp 


iu 

2 
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j  «'  J.  (^)  ^=^ 


-iuO 


2a 


dO  . 
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For  small  x,  J  (x)^^x     and  we  see  that  I,  and  I^  are  of  lower  order  in 
'  n  12 

a  than  I  ,  so  that  these  integrals 
o 

Hence  (2-50)  and  (2-31)  reduce  to 


a  than  I  ,  so  that  these  integrals  may  be  neglected  in  comparison  with  I  . 


(^^.2) 


e  =  h 
X    y 


LAI  , 
o  ' 


{^■3) 


e  =e  =h  =h  =0 
y    z    X    z 


Thus  in  a  system  with  a  small  angiilar  aperture,   the  field  in  the  image 

region  is  effectively  linearly  polarized;  and  the  directions  of  the  two 

field  vectors  in  the  image  region  axe  the  same  as  their  directions  in 

the  object  space.  We  see  that  in  this  case  the  image  field  is  conrpletely 

specified  by  one  (complex)  scalar  wave  (e  or  h  ),  which  is  represented  by 

X    y 

the  integral  I  alone. 

To  evaluate  I  we  introduce  a  new  vairLable  p  =  ©/a  and  obtain 
o 


(k.k) 


1   (u,v)  =  aa  exp 


iu 


a 

•-    -^   0 


j  p  J^  (vp)  exp 


1.   2 
2IUP 


dp 


The  integral  on  the  right  may  be  evaluated  in  terms  of  two  of  the  Lommel 
functions 


(^.5) 


cp     _    n  +  2s 

U  (u,v)  =T~   (-1)®  {^)  J  ^  ^  (v)  , 

n       ^—7:  V        n  +  2s 

8=0 


introduced  by  Lommel  in  the  aneilysis  already  referred  to. 
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In  fact* 


J. 
ih.6)  j      p  J^(vp)  exp     -  2  iup       dp  =     -  exp     -  ^  iu 


U^(u,v)   +   iU2(u,v)j 


From  (^.2),    (U.U)   and  (^.6)  we  obtain  the  following  expressions  for  the  only 
two  non-vanishing  field  conrponents  in  the  image  region 


(^.7) 


exp 


a 


f 


e     =  h 
X         y 


2iAa 


U^(u,v)   +  iU2(u,v) 


u 


The  expressions  (5.1^)  for  the  energy  densities  now  beccme 


(h.Q)  <w>  +  <w  > 

^  e     ^  m 


|<w> 


U^^(u,v)  +  u|(u,v)J 


*See  Lommel  (1885)  or  Bom  and  Wolf  (1959),  S  8.8.  The  U^-series  converge 

for  «n  u  and  v  values,  but  are  convenient  for  computations  only  when 

|u/v|  <  1.  When  |u/v|  >  1,  (h-.k)   may  be  evaluated  in  terms  of  two  of  the 

Lomnel's  V  -functions  defined  by 
n 


00 


V„(u,v)  =  Yl     (-1) 

"  8=0 


/  \  n+2s 

(0 


Vas  (''  • 


The  U  and  V-functions  are  related  by  means  of  the  formula   [cf .  Watson, 
1952,   557] 


vn-1 


"n+1  ^'''''^  ■   (-1^         Vl  ^"^'""^  =  ^^°    ^ 


1  v2 

^u-H-.n* 
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The  formulas  (3«2l)  for  the  components  of  the  Poynting  vector 


beccme 


(".9)     <S^>  =  (B^)    =  0 


<=.>-  -4ki''   #  -.    [u/Kv).u|(.,v)]. 


Thus,  vhen  a  Is  small,  the  energy  flov,  as  represented  by  the  time  averaged 
Poynting  vector,  is  in  the  direction  of  the  positive  z-axLs.  On  ccanparison 
of  C+.S)  and  ('♦■.9)  we  also  see  that  the  relation 


(i^.lO)     I  <S)  I   =  c  (v) 


then  holds  everywhere  in  the  image  region. 

The  expressions  on  the  right  of  (^.8)  and  (^.9)  are  proportional 
to  the  cleissical  solution  of  Loramel  and  Struve  for  the  'intensity'  in  the 
image  region  of  an  aberration- free  system. 

When  u  =  0  (focal  plane)  we  have  from  (^.7)  since 


U^(u,v)     J^(v)  U2(u,v) 

—  =  ,.   Lim 

u  V    ' 

u->  0  u-^  0 


Lim   —  =  .   Lim   =  0  , 

U  V     '  u 


2  /2Ji(v)^ 
(1^.11)    e^  =  hy  =  -  iAa''  ' 


and  i^.B)   and  (^.9)  become 


(^•^)  <-e>    =    <-»)   =  I  <')    -    55  I  (=)  I  =    #  (^ 
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The  expression  on  the  right  of  (^.12)  is  proportional  to  the  classical 
solution  of  Airy  for  the  'intensity'  in  the  image  plane  of  an  aberration- 
free  system. 

5.   Images  formed  by  an  unpolarlzed  wave 

So  far  we  have  assumed  that  the  vave  entering  the  system  is  linear- 
ly polarized.  We  shall  now  briefly  consider  the  case  when  the  wave  is  un- 
polarized. 

An  unpolarized  incident  wave  may  beccane  partially  polarized  on 
refraction  (or  on  reflection)  at  the  successive  surfaces  of  the  system. 
However,  if  the  angles  of  incidence  at  each  surface  are  small  as  is  aasumed, 
the  degree  of  polarization  introduced  in  this  way  will  also  be  small  and  we 
may,  therefore,  consider  the  wave  emerging  from  the  system  to  be  effectively 
unpolarized.  If  further,  we  assume  that  the  wave  is  quasi -monochromatic, 
with  mean  angular  frequency  cq,  the  expressions  for  the  time  averaged  energy 
densities  and  energy  flow  may  be  obtained  by  averaging  the  corresponding 
expressions  relating  to  the  polarized  wave  over  all  possible  states  of  pol- 
eoT-zation. 

To  carry  out  this  averaging  it  is  convenient  to  choobe  a  new  set 
of  coordinate  axes  0|,  Ot]  at  right  angles  to  each  other  and  to  the  axis  of 
revolution.  Let  0  be  the  angle  between  0|  and  the  direction  Ox  of  the 
electric  vibrations  of  the  incident  wave,  and  let  0p  be  the  angle  between 
0|  and  the  azimuth  which  contains  the  point  P  of  observation  (see  Figure  6). 

Then  the  first  formula  in  (5.1^)  for  the  time  average  electric 
energy  density  at  a  point  P  of  observation  beccsnes,  since  0  =  0p  -  0^  , 

(5-1)    <-e)  =1^   llol'-^|lll'<=o«'(V«^o)^l^2l'^2cos[2(0p.0j]^(V2*) 

If  we  denote  by  verticsJ.  bar  an  average  over  alT  possible  values  of 

0  (0  <  0  <  2jt),  it  follows  that  when  the  incident  wave  is  unpolarized 
''o^  -  ^o 
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FIGURE  6 


For  caption  to  Figure  see  page  i^L« 
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the  time  averaged  electric  energy  density  at  P  is 

We  note  that  /v  \  is  independent  of  0p  so  that  the  loci  of  the  constant 
averaged  electric  energy  are  surfaces  of  revolution  about  the  axis  of  the 
system,  as  vas  to  be  expected  from  symmetry.  In  a  strictly  similar  manner, 
we  obtedn  expressions  for  /w  \  and  /wS   ,  and  we  have  in  all 

(5-5)     (?)  '    K)    '-2(^}    '4    {  l^cl'  *  ^l^ll'  *  ^h\"  I  • 

Comparison  with  the  last  expression  in  {j>.l^)   shows  that  each  term  in  (5'5) 
is  also  equal  to  one  half  of  the  total  energy  density  of  the  polarized  wave. 
The  distribution  of  this  quantity  in  the  focal  plane  has  already  been  given 
in  Figure  5. 

Next  consider  the  energy  flow.  It  is  agEiin  convenient  to  trans- 
form first  the  expressions  for  the  polarized  case  so  that  the  coordinates 
which  appeeu:  in  them  are  referred  to  the  fixed  axes  Oi,   Ori,  independent  of 
the  direction  of  polarization  of  the  incident  wave.  Let  /s,\  and  {s  ) 
be  the  components  of  the  time  averaged  Poynting  vector  in  the  directions 
05,  Ot).  Then 

(5  ,^)     .  ih)    -    <^x>  -«  ^o  -  <^y>  «^  ^o  ' 
<0  =  <Sx)  «^°^o^  <0  =°«^o  • 


*  The  factor  £    entering  the  constant  A  =  rtf^  /X  may  now  be  related  to  the 

O  r r         O 

average  electric  energy  density  (v  ^  .      of  the  incident  wave  in  the  object 
space  by  means  of  the  formula 


<'e)i-  is  (?o-V>   =  m:  K 


3$ 


On  substituting  for  -C  S  >  and  <[  S  ">  from  (5. 21),  where  again  0  is 

X        y 

replaced  by  0  -  0  ,  ve  find  that  ^S.y    and  *Cs  >  are  given  by  the  same 

fonmilas  as  <  S  >  and  <.  S  >  but  vlth  0^  written  in  place  of  0.  Since 

X  ^        J 

these  formulas  are  independent  of  0  ,  they  apply  not  only  to  the  polar- 
ized wave  but  also  to  the  unpolarized  one,  i.e. 


(5.5a) 


(5.5b) 


The  remaining  component,  in  the  direction  of  the  axis  of  revolution  iSj 
as  before,  given  by  the  last  formulas  in  (5-21)  , 


(5.5c)  <\>  =<s^>  =s- il^cl  -1^2!' 


It  immediately  follows  that  the  resiilts  expressed  by  (3-22),  (5'23)  and 
(5.24)  hold  also  for  the  unpolarized  wave. 

When  the  angular  semi-aperture  a  is  sufficiently  small  (a  — ^O)  one  may, 
as  in  ^  4,  neglect  I^  and  I^  in  comparison  with  I  .  The  above  formulas  then 
become,  when  (h-.k)   and  (4.6)  are  used: 

(5.6)     ^>  =  T^y    =    ^<^   =  ^  ^  [u/(u,v)  +  u/(u,vT]  , 


(5.7a)    <s^>  =  <S^>=  0, 

.2  k 


(5.7b)    <S^>  =  ^  \   [u/(u,v)  +  u/(u,v7]  . 


Here  U  and  U  are  again  two  of  the  fimctions  defined  by  (4.5), 
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For  points  of  observation  in  the  focal  plane    (u  =  O)   these  formulas 
r«duce  to 


2  U 


(5.8)  <v^>    =<w^>   =|<v>  =  2^  i<f>|   -^ 


As  already  noted  in  connection  with  the  corresponding  formulas  (k.Q),    (k.^) 
and  {k.l2),   the  quantities  on  the  right  hand  sides  of  (5'6)  sni   (5.7b)  are 
proportional  to  the  classical  solution  of  Lommel  and  Struve  for  the  'inten- 
sity' in  the  image  region  of  an  aberration- free  system,  and  the  expression 
on  the  right  hand  side  of  (5'8)  is  proportional  to  the  classical  solution 
of  Airy  for  the  'intensity'  in  the  image  plane. 

Let  us  denote  by  suffix  zero  expressions  (such  as   (5-6))  which  refer 
to  a  low  aperture  system  (a  -^  O)  .  Then  from  (^.8),  (^.9), (5-5),  (5'6) 
and  (5'7)  the  following  relations  are  seen  to  hold: 


(5-9)  <S^>o  =  <S,>o  =  '=<^>o<=<^>  =  '=<^> 


(5.10)  <w^>^  =   <w^>^  =  |<w>^  <  |<v>    =  <v^>   =  <v^>       . 

In  piractice,  detectors  of  electric  energy  are  usually  used.  The  relation 
<C^  y      —  ^^  ^  implies  that  the  pattern  recorded  in  any  particuleir 
receiving  plane  (u  =  const)  will  then  be  broader  in  a  system  with  a  wide 
angiilar  aperture  than  in  one  with  a  low  angular  aperture.  For  patterns  in 
the  focal  plane,  this  effect  is  seen  in  Figure  5  which  also  shows  that 
<^w  N  has  no  exact  zeros,  nor  pronounced  subsidiary  maxima. 

Finally,  some  of  the  main  data  which  relate  to  the  structure  of  the 
image  in  the  focal  plane  of  an  aplanatic  system  are  summcarized  in  a  Table 
below.  Data  relating  to  both  polarized  and  unpolarized  incident  waves  are 
given . 
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TABLE 


Value 

Quantity 

a-*0 

a  =50° 

a  =60° 

V 

a.  if 

a -^90° 

<^v^(0,v,0)> 

=    <-J°'-'|)> 

0.93 

0.98 

1.16 

1-55 

1.1+8 

80 

<w^(O,v,0)> 

=    <wj0,v,0)> 
=    <wJO,v,0)> 

0.95 
0.93 

0.95 
0.95 

0.89 
1.00 

0.86 
1.05 

0.81+ 
1.05 

<;s^(O,v,0)> 

=  <S^(O,v,0)> 

0.93 

0.95 

0.93 

0.93 

0.91 

^v^(0,v,0)> 

=    <wJO,v,|)> 

l.J+0 

iM 

1.73 

1.95 

2.10 

60 

<w^(0,v,|)> 
<yg(O,v,0)> 

=    <wJO,v,0)> 

i.Uo 
l.Vo 

1.58 
1.1+2 

1-52 
1.50 

1.30 
1.55 

1.26 

=    <Wj^(O,v,0)> 

1.58 

<S^(O,v,0)> 

=    <Sj,(O,v,0)> 

i.40 

1.1+0 

1.39 

1.58 

1.56 

<_v^(0,v,0)> 

=   <vJO,v4)> 

1.85 

1.9^ 

2.27 

2.50 

2.65 

ko 

<w,(0,v,|)> 
^w^(O,v,0)> 

=    <wj0,v,0)> 

1.85 
1.85 

1.82 
1.87 

1.70 
1.99 

1.70 
2.05 

1.61+ 

=    <wJO,v,0)> 

2.10 

<S^(O,v,0)> 

=  <S^(O,v,0)> 

1.35 

1.83 

1.82 

1.82 

1.79 

^w^(0,v,0)> 

-    <-n.(°'^'|)> 

2.50 

2.50 

2.93 

5.1'^ 

5.25 

20 

^w^(O,v,0)> 

=    <v  (0,v,0)> 
=    <wJO,v,0)> 

2.56 
2.36 

2.55 
2.1+1 

2.22 
2.56 

2.16 
2.66 

2.10 
2.75 

<^S^(O,v,0)> 

=    <S^(O,v,0)> 

2.56 

2.56 

2.3'+ 

2.53 

2.27 

^w^(0,v,0)> 

=    <-J°'-'|)> 

3.85 

1+.15 

I+.75 

1+.75 

I+.70 

l£t 

<-e(°'^'l^> 

=    <vJO,v,0)> 

5.83 

5.77 

3.56 

5.1+0 

5.29 

Min. 

<w^(O,v,0)> 

=    <wJO,v,0)> 

5.83 

5-95 

- 

- 

- 

<S^(O,v,0)> 

=    <S^(O,v,0)> 

3.85 

5.80 

5.80 

5.76 

5.70 

<^w^(0,v,0)> 

=    <-J°^-'|)> 

7.02 

7  05 

8.10 

8.10 

8.00 

■  and 

<^v^(O,v,0)> 

=    <wJO,v,0)> 

7.02 
7.02 

6.95 

6.85 

6.70 

6.55 

Min. 

=    <wJO,v,0)> 

- 

<S^(O,v,0)> 

=    <S^(O,v,0)) 

7.02 

7.02 

7.02 

7.02 

6.95 

Caption  to  Table   is  on  page  1^2. 
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eAPTIONS  TO  FIGURES 


Figure  1  :     The  meridional  plane  of  a  ray.  The  axis  Ox  is  in  the 

direction  of  the  electric  vector  e  in  the  object  space. 

Figure  2  :     (a)  Contours  of %(v,0),  (b)  Contours  of  p(v,0).  Polar- 
ization of  the  electric  field  in  the  focal  plane  of  an 
aplanatic  system  of  angxolar  semi-aperture  a  =  60°  on 
the  image  side. 

The  field  in  the  object  space  is  linearly  polarized 
with  its  electiric  vector  in  the  azimuth  0=0.  The  field 
in  the  focal  plane  is,  in  general,  ellipticaJLly  polarized, 
with  one  axis  (q)  in  the  focal  plane,  and  the  other  axis 
(p)  pe]T)endicular  to  the  focal  plane.  ^(v,0)  is  the 
angle  between  the  plane  of  the  polarization  ellipse  of 
the  electric  vector  and  the  meridional  pleuie  0=0,  emd 
p(v,0)  is  the  ratio  of  the  p  and  q  axes.  Radial  distances 
are  measured  in  v-units.  The  dotted  lines  indicate  the 
dark  rings  of  the  Airy  pattern  represented  by  ('+.12). 

Figiire  3  :     Contours  of  the  time  averaged  electric  energy  density 
<^w  y     in  the  focal  plane  of  an  aplanatic  system  of 
angular  semi-aperture  a  on  the  image  side. 

The  field  in  the  object  space  is  linearly  polarized 
with  its  electric  vector  in  the  azimuth  0=0.  The  crosses, 
dots  and  circles  indicate  maxima,  minima  Eind  saddle  points 
respectively.  Radisil  distances  are  measured  in  v-units. 
The  values  are  normalized  to  100  at  the  focus.  Figure  (a), 
which  represents  the  limiting  case  a  — >  0  is  identical 
with  the  classical  intensity  pattern  of  Airy,  characterized 
by  (J+.12).   The  dotted  circles  in  Figures  (b)  -  (e) 
represent  the  dark  rings  in  the  Airy  pattern. 
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Figure  k-  :  (a)  Meridional  section  0  =  0.   (b)  Meridional  section  0  =  ^  . 

The  variation  of  the  time  averaged  electric  energy    /v  \ 
along  the  two  principal  meridional  section  of  the  focal  plauae 
in  an  aplanatic  system  of  selected  angtilar  semi -aperture  a  in 
the  image  space. 

The  field  in  the  object  space  is  linearly  polarized  with 
its  electric  vector  in  the  azimuth  0=0.  The  values  are 
normalized  to  100  at  the  focus.  The  curves  representing  the 
limiting  case  a  ->  0  are  identical  with  the  classical  Airy 
intensity  curve  given  by  (^.12). 

Figure  5  :     The  variation  of  the  time  averaged  total  energy  density 
(Figure  a)  and  of  the  only  non- vanishing  ccnrponent  of  the 
time  averaged  energy  flow  (Figure  b)  along  any  meridional 
section  0  =  const,  of  the  focal  plane  of  an  aplBnatic  system 
of  angular  semi-aperture  a  on  the  image  side.  The  curves 
represent  the  case  of  polarized  as  well  as  unpolarized  wave. 
The  cujrves  in  (a)  also  represent  the  variation  of  the  time 
averaged  electric  and  magnetic  energy  densities  of  an 
unpolarized  wave.  The  values  are  normalized  to  100  at  focus. 

Figure  6  :     Change  of  reference  system.  Section  at  right  angles  to  the 
axis  of  revolution. 
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Caption  to  Table 


Comparison  of  data  relating  to  the  structure  of  the  electro- 
magnetic field  in  the  focal  plane  (u  =  O)  of  an  aplanatlc  system  of 
selected  angular  semi-aperture  a  on  the  image  side. 

<^w  (u,v,0)^  and  <^w  (u,v,0)^  represent  the  time  averaged 
Electric  energy  density  and  magnetic  energy  density  respectively. 
<Cs  y    represents  the  only  non- vanishing  component  (in  the  direction  of 
the  axis  of  the  system)  of  the  time  averaged  energy  flow.  All  the  quan- 
tities are  normalized  to  100  at  the  focus  u  =  v  =  0. 

Quajitities  with  a  vertical  bar  refer  to  an  unpolarlzed  wave, 
the  others  to  linearly  polarized  wave  in  the  object  space,  with  its 
electric  vector  In  the  azimuth  0  =  0.   The  tlilrd  row  in  each  division  of 
the  table  also  represents  one  half  of  the  time  averaged  total  energy 
density  <l,w  ">  ,  i.e. 


<-e>  =  <"m>  =!<->  =  !<->  • 


The  entries  in  the  column  a — >  0  are  identical  with  values  given  by  the 
intensity  formulas  (4.12)  of  Airy. 
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3700  East  Pont  lac  Street 
Fort  Wayne  1,  Indiana 
Attn:  Technical  Library 

Federal  Telecommunication  Labs. 
?00  Washington  Avenue 
Nut ley  10,  New  Jersey 
Attn:  Technical  Library 

The  Gabriel  Electronics 
Division  of  the  Gabriel  Company 
13?  Crescent  Road 
Needham  Heights  9U,  Mass. 
Attn:  Mr.  Steven  Galagan 


General  Electric  Advanced  Electronics  Center 
Corr.ell  University 
Ithaca,  New  York 
Attn:   J.  B.  Travis 

General  Electric  Company 

Elpctronics  Park 

Syracuse,  New  York 

Attn:   Documents  Library,  B.  Fletcher 

Building  3-'lb3A 

General    Precision  Laboratory,    Inc. 

63  Bedford  Road 

Pleasantvllle,    New  York 

Attn:      Mrs.   Mary  G.   Herbst,    Librarian 

Goodvear  Aircraft  Corp. 

1210  Masslllon  Road 

Akron  15,   Ohio 

Attn:     Library  DA20  Plant  A 


Granger  Associates 
Electronic  Systems 
966  Commercial  Street 
Palo  Alto,    California 
Attni      John   V.   N.   Granger, 


President 


Grumman  Aircraft  Engineering  Corporation 
Bethpage,  Long  Island,  N.  Y. 
Attn:  Mrs.  A.  M.  Gray,  Librarian 

Engineering  Llbrar-/,  Plant  No.  5 

The  Halllcrafters  Company 

UbOl  West  5th  Avenue 

Chicago  2b,  Illinois 

Attn:  La Verne  LaGloia,  Librarian 

Hoffman  Laboratories,  Inc. 
3761  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Microwave  Laboratory 
Building  12,  Room  2617 
Culver  City,  California 
Attn:  M.  D.  Adcock 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 

Attn:  Dr.  L.C.  Van  Atta,  Associate  Director 
Research  Labs. 

Hycon  Eastern,  Inc. 
75  Cambridge  Parkway 
Cambridge,  Mass. 
Attn:  Mrs.  Lois  Seulowltz 
Technical  Librarian 

International  Business  Machines  Corp. 

Military  Products  Division 

590  Madison  Avenue 

New  York  33,  New  York 

Attn:  Mr.  C.F.  McElwsin,  General  Manager 

International  Business  Machines  Corp. 
Military  Products  Division 
Owego,  New  York 

Attn:  Mr.  D.  I.  Marr,  Librarian 
Department  b59 

International  Resistance  Company 
bOl  N.  Broad  Street 
Philadelphia  8,  Pa. 
Attn:   Research  Library 

Jansky  and  Bailey,  Inc. 
1339  Wisconsin  Avenue,  N.  W, 
Washington  7,  D.  C. 
Attni  Mr.  Delmer  C.  Ports 

Dr.  Henry  Jasik,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  Hew  York 

Electromagnetic  Research  Corporation 
711  Ibth  StJ'eet,  N.  W, 
Washington  5,  D.  C. 
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Lockheed  Aircraft  Corporation 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-75,  Plant  A-1,  Bldg.  63-1 

Purbank,  California 

Attn!   N.  C.  Harnols 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn!  Mr.  Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attn  I  Engineering  Library 
Antenna  Design  Group 

Maryland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attn!  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 

The  W.  L.  Maxscn  Corporation 

Ii60  West  ^Uth  Street 

New  York,  N.  T. 

Attn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,    St.   Louis   3,   Missouri 

Attn:      R.   D.   Detrich,    Engineering  Library 

McMillan  Laboratory,    Inc. 

Prownville  Avenue 

Ipswich,    Massachusetts 

Attn:     Security  Officer.    Docuirent  Room 

Melpar,    Inc.  (2' 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

Attn!     Engineering  Technical  Library 

Microwave  Development  Laboratory 
90  Broad  Street 
Babson  Park  57,    Massachusetts 
Attn:     N.   Tucker,    Oeneral  Manager 

Microwave  Ralietion  Comoany  Inc. 

19223  South  Hamilton  Street 

Gardena,    California 

Attn:     Mr.   Morris   J.  Ehrlich,    President 

Chance  Vought  Aircraft,    Inc. 
931I4  West   Jefferson  Street 
Dallas,   Texas 
Attn:      Mr.   H.   3.   White,    Librarian 

Northrop  Aircraft,    Inc. 
Hawthorne,    California 

Attn:     Mr.  E.   A.   Freitas,    Library  DeptJlli5 
1001  E.  Broadway 

Remington  Rand  Univ.   -  Division  of  Sperrv 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,    Pennsylvania 
Attn:     Mr.    John  F.  McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,    Inc. 
12211!  Lakewood  Boulevard 
Downey,   California 
Attn:     Engineering  Library  U95-115 

North  American  Aviation,    Inc. 
l^s  Angeles   International  Airport 
Los  Angeles  Ii5,    California 
Attn:     Engineering  Technical  File 

Page  Commi!nic9tions  Engineers,    Inc. 
710  Fourteenth  Street,    Northwest 
Washington  5,    D.   C. 
Attn:     Librarian 

Philco  Corporation  Research  Division 

Branch  Library 

U700  Wissachickon  Avenue 

Philadelphia  hh.   Pa. 

Attn:     Mrs.  Dorothy  S.  Collins 


Plckard  and  Burns,    Inc. 

2bO  Highland  Avenue 

Need  ham  9lj,   Mass. 

Attn:     Dr.    J.  T.   DeBettencourt 

Polytechnic  Research  and  Development 

Company,  Inc. 

202  Tlilary  Street 

Brooklyn  1,  New  York 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 
Attn:  Dr.  John  Ruze 

Radiation,  Inc. 

P.  0.  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 


Radio  Corp.  of  America 
RCA  Laboratories 
Rocky  Point,  New  York 
Attn:   P.  S.  Carter,  Lab. 


Library 


RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,  New  Jersey 
Attn!  Miss  Fern  Cloak,  Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Prcmcts 
Building  10,  Floor  7 
Camden  2,  New  Jersey 
Attn:  Mr.  Harold  J.  Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Adm"!  nlstrator 

The  Ram^'-Wooldridge  Corporation 
P.'^.   ''ox  Ij5Ii53   Alrnort   Station 
Los  Angeles  )i5,    California 
Attn:     Margaret  C.   Whitnah, 
Chief  Librarian 

Hoover  Micowave  Co. 
9592  Baltimore  Avenue 
College  Park,    Maryland 

Director,    USAF  Project  RAND 

Via:     Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  "ain  Street 

Santa  Monica,    California 

Rantec  Corporation 

Calabasas,   California 

Attn:      Grace  Keener,   Office  Manager 

Raytheon  Manufacturing  Company 
Kissila  Systems  Division 
Bedford,    Macs. 
Attn:     Mr,    Irving  Goldstein 

Raytheon  Manufacturing  Company 
Wayland  Laboratory,    State  Road 
Way land,    Mass. 
Attn:     Mr.   Robert  Borts 

Raytheon  Manu'"acturing  Comrany 
Wayland  Laboratory 
Wayland,   Mass. 

Attn:     Miss  Alice  0.  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,    Long  Island,   N.  Y. 
Attn:     Engineering  Library 

Thru:     Air  Force  "lant  Representative 
Republic  Aviation  Corp. 
Farmingdale,   Long  Island,   N.Y. 

Rheen  Hanu*'acturing  Conpany 
9236  East  Hall  Road 
Downey,    California 
Attn:      J.  C.    Joerger 

Trans-Tech,    Inc. 
P.  0.  Box  3U6 
Frederick,   Maryland 


Ryan  Aeronautical  Comoany 
Lindbergh  Field 
San  Diego  12,    California 
Attn:     Library  -  unclassified 

3age  Laboratories 
159  Linden  Street 
V/ellesley  fil,   Mass. 

.Sanders  Associates 
95  Canal  Street 
Nashua,   New  Hampshire 
Attn:     N.   R.  Wild,    Library 

Sandia  Corporation,   Sandia  Base 

P.O.  Box  5''00,    Albuquerque,    New  Mexico 

Attn:     Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,    Long   Island,    New  York 

Attn:      Florence  W.  Turnbull,    Engr.   Librarian 

Stanford  Research  Institute 

Menlo  Park,    California 

Attn:     Library,   Engineering  Division 

Sylvania  Electric  Products,    Inc. 
100  first  Avenue 
Walthai:  51i,    Mass. 

Ittn:     Charles  A.  Thornhill,    Report  Librarian 
Walthan  Laboratories  Library 

Systems  Laboratories  Corporation 
U4P52   Ventura  Boulevard 
Sherman  Oaks,    California 
Attn:     Donald  L.   Margerum 

THG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y. 

Attn:  M.  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

161  Devonshire  Street 

Boston  10,  Mass. 

Attn:   A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Chu  Associates 
P.  0.  Box  3.''7 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Company 
Missil"  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  Laboratory 
225  Crescent  Street 
Waltham,  Mass. 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50b,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Cnief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Dr.  W.  C.  Hoffman 

Commander 

AF  Office  of  Scientific  Research 

Air  Research  and  Development  Command 

lUth  Street  and  Constitution  Avenue 

Washington,  D.  C. 

Attn:  Mr.  Ottlng,  SRY 

Westlnghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int '1  Airport  Box  7li6 
Baltimore  3,  Maryland 
Attn!  Engineering  Library 
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Wheeler  Laborstorifts,    Inc. 
122  Cutter  Mill  Road 
Oreat   Meek,    New  York 
Attn:     Mr.   Harold  A.   Wheeler 

Zenith  Plastica  Co. 
Box  91 

fiardena,   California 
Attn:     Mr.   S.   S.   Oleestcy 

Library  Geophysical   Institute 
of  the  University  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  li,   California 

Attn:     Dr.  Sainuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elee.  Eng.   Electronics 
Research  Lab. 

University  of  California 
Electronics  Research  Lab. 
332  Corv  Hall 
Berkeley  h,   California 
Attn:     J.   R.  Whlnnery 

California   Institute  of  Technology 
Jet  Propulsion  Laboratory 
IjPOO  Oak  Grove  Drive 
Pasadena,   California 
Attn:     Mr.    T.   E.   Hewlan 

California  Institute  of  Technology 
1201  E.   Calif  omi  a  Street 
Pasadena,   California 
Attn:     Dr.  C.    Papas 

Carnegie  Institute  of  Technology, 
Schenlev  Park 

Pittsburgh  13,    Pennsylvania 
Attn:     Prof.   A.   E.   Helns 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,   New  York 

Attn:     Prof.  G.   C.   Dalman 

Tniverslty  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,    Florida 

Attn:      Prof.   K.   H.   Latour,    Library 

Llhr-ary 

Georgia   Institute  of  Technology 

Engineering  Experiment   Station 

Atlanta,    Georgia 

Attn:     Mrs.    .'.H.   Crosland,    Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,    303A   Pierce  Hall 

Oxford  Street,   Cambridge  3'',    Mass. 

Attn:     Mrs.   E.L.   Hufschmidt,   Ubrarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,   Mass. 
Attn:     Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,    Illinois 

University  of  Illinois 
College  of  Enginearlng 
Urbana,    Illinois 

Attn:      Dr.    P.   E.   Moyes,    Department  of 
Electrical    Engineering 

The  Johns  Hopkins  University 
Homewood  Campus 
Department  of   Physics 
Baltimore  IQ,   Maryland 
Attn:     Dr.   Donald  E.  Kerr 

Sandla  Corporation 
Attn:     Organization  l'i23 
Sandla  Base 
AlbuouerfTue,   New  Mexico 


Applied  Physics   Laboratory 
The  JohrS Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,   Maryland 
Attn:     Mr.    George  L.   Seialstad 

Massachusetts   Institute  of  Technologv 

Research  Laboratory  of  Electronics 

Room  20B-221 

Cambridge  39,  Massachusetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
Electronic  Defense  Croup 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J.  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

Universitv  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilanti,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  lli,  Minnesota 

Attn:  Mr.  Robert  H.  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:   R.  E.  Beam 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tlce 

Dept.  of  Elec.  Engineering 

The  UnlTersity  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attn:  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A.  Ollner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
S5  Johnson  Street 
Brooklyn,  New  York 
Attn;  Mr.  A.  E.  Laammel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  o.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:  Prof.  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Hinglns 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G.  Held,  Associate  Professor 

Stanford  University 
Stanford,  California 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agricultxire 

and  Mechanic  Arts 

State  College,    New  Mexico 

Attn:     Mr.   H.  W.   Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn;  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof.  S.  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momingside  Heights 

New  York,  N.  Y. 

Attn:  Dr.  Schlesinger 

McGill  University 

Montreal,  Canada 

Attn:  Prof.  Q.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schultj 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsvlvanla 

Institute  of  Cooperative  Research 

3liOO  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept.  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 
W.  Cumberland  Avenue 
Knoxvllle  16,  Tennessee 

Universitv  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  tinlversity 

Detroit,  Michigan 

Attn:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C.  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parti 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  U,    Pennsylvania 


Ionosphere  Research  Laooratory 

Pennsylvania  State  College 

State  College,  Penneylvanla 

fT"!:  Professor  A.  H.  Waynick,  Director 

Institute  of  Mathematical  Sciences 

25  Waverly  Place 

New  York  3,  Ncvi  York 

ATTNt  Librarian  (3) 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
ATTNi  Dr.  Robert  Kalaba 

National  Bureau  of  Standards 

Washington,  D.  C, 

ATTN:  Dr.  W.  K.  Saunders 

Applied  Mathematics  and  Stati.sti'^s  Lab, 

Stanford  Oniverslty 

Stanford,  (ill<'om'ia 

ATTN:  Dr.  Albert  H.  Bowker 

Departme'it  of  Physics  and  Astrjnomy 
Michigan  State  College 
East  Lansing,  Michigan 
ATTM:  Dr.  A.  Leltnor 

University'  of  Tennessee 
Knoxville,  Tennessee 
ATTN:  Dr.  Fre  I  A.  Ficken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN:  Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Diego  12  California 

Attn:   Dr.  Edward  Oerjuoy 

Department  of  Physics 
Amherst  College 
Amherst,  Mass. 
ATTN:  Dr.  Arnold  Arons 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
ATTN:  Dr.  A.  Erdelyi 

Mathematics  Derertment 

Stanford  nnl versity 
Stanford,  California 
ATTN:  Dr.  Harold  Levlne 

University  of  Minnesota 

Minnpapolis  I'l,  Minnesota 

ATTN:  Professor  Paul  C.  Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,  California 

ATTN:  Professor  Bernard  Epstein 

Aoplied  Physics  Laboratory 
The  JohnP  Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring,  Mar:;-land 
ATTN:  Dr.  B.  S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,  D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technolo^ 
Cambridge  39,  Mass. 
ATTN:  Dr.  L.  J.'  Chu 

Nuclear  Development  Associates,  Inc. 

5  New  Street 

'Aite  Plains,    New  York 

ATTN:  Library 

California  Institute  of  Technology 

Electrical  Engineering 

Pasadena,  California 

ATTN:  Dr.  Zohrab  A.  Kaprielian 


Dr.  Rodman  Doll 
311  W.  Cross  Street 
Ypsilanti,  Michigan 

California  Inst,  of  Technology 
Pasadena,  California 
ATTN:   Mr.  Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Mass. 

National  Bureau  of  Standards 
Boulder,  Colorado 
ATTN:  Dr.  R.  Gallet 

Dr.  Solomon  L.  Schwebel 

3689  Louis  Road 

Palo  Alto,  California 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  lb,  Minnesota 
ATTN:  Exchange  Division 

Department  of  Mathematics 
University  of  California 
Berkeley,  California 
ATTN:  Professor  Bernard  Friedman 

Lincoln  laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Bex  73 

Lexington  73,    Massachusetts 

ATTN:     Dr.   Shou  Chin  Wang,    Room  C-351 

Melpar,    Inc., 

3000  Arlington  Boulevard 

Falls  Church,    Virginia 

ATTN:     Mr.   K.   S.   Kelleher,    Section  Head 

Hq.   Air  Force  Cambridge  Research  Center 

Laurence  G.    Hanscom  Field 

Bedford,    Mass. 

ATT-J:     Mr.   Francis  J.   Zucker,   CSRD 

Hq.   Air  Force  Cambridge  Research  Center 

Laurence  G.   Hanscom  Field 

Bedford,    Mass. 

ATTN;     Dr.   Philip  Newman,    CRRK 


Smyth  Research  Associates 
3555  Aero  Court 
San  Diego  3,   California 
ATTN:     Dr.    John  B.  Smyth 

Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,    California 
ATTN:     Dr^orgea  G,  Weill 

Naval  Research  Laboratory 

Washington  25,    D.  C. 

ATTN:      Henry  J.   Passerini,   Code  5278A 

Dr.   George  Kear 
5  Culver  Court 
Orinda,   California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,    New  York 

ATTN:     Dr.  Nathan  Marcuvltz 


Mr.  N.  C.  Gerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Barrar 
Systems  Development  Corp. 
21:00  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

II45  Palisade  Street,  Dobbs  Ferrv,  N.  Y. 
ATTN:  Dr.  N.  H.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland 
ATTN:  Dr.  Elliott  Montroll 

Depai-tment  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5,  Mo. 

ATTN:   Professor  J.  Van  Bladel 

Department  of  the  Navy 

Office  of  Naval  Research  Branch  Office 
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